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Abstract 

The ion trap quantum computer proposed by Cirac and Zoller [Phys. Rev. Lett. 74, 
4091 (1995)] is analyzed for decoherence due to vibrations of the ions. An adiabatic 
approximation exploiting the vast difference between the frequencies of the optical 
intraionic transition and the vibrational modes is used to find the decoherence time 
at any temperature T . The scaling of this decoherence time with the number of ions 
is discussed, and compared with that due to spontaneous emission. 



I. INTRODUCTION 

Since Shor's discovery [1] of an algorithm for factorization of a composite number 
of order 2 L in ~ L 3 steps on an ideal quantum computer (QC), a great deal of 
effort has gone into developing quantum computational theory, and related ideas 
in information transfer and cryptography. This work has shed new light on both 
quantum mechanics and computational complexity theory. As Landauer [2] has so 
pungently said, however, writing down a Hamiltonian is not the same as specifying an 
apparatus, and it is also necessary, at some time, to look into building real machines. 
Landauer's criticism has not gone completely unanswered, and concrete proposals for 
implementing QC's have been put forward. Perhaps the most promising of these is 
by Cirac and Zoller (CZ) [3] - but see also Ref. [4]. It is based on a linear array of 
trapped ions driven by a precisely timed sequence of laser pulses. This proposal is 
being taken seriously enough that at least one group in the world is trying to build a 
prototype [5]. 

The CZ QC consists of N identical ions, trapped and cooled in a linear rf Paul trap. 
The ions form a linear array with nonuniform spacings determined by their mutual 
Coulomb repulsion and the effective trapping potential (see Fig. 1). Two internal 
states of each ion, \g) and |e), serve as the quantum bit, and laser pulses {tt, ir/2, 
etc.) drive g «-> e transitions and thus implement one-bit gates. These transitions 
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can be driven either by lasers tuned to the direct g <-> e transition frequency u>q, or 
by Raman pulses in a A system, where the lower levels form the quantum bit, and are 
chosen from the ground multiplet of the ion to minimize spontaneous emission. The 
innovative idea due to CZ is to use the center of mass longitudinal vibrational mode 
of the array as a bus that enables the execution of two-bit gates. This can be done 
by a combination of one-ion pulses and pulses detuned by u) z , the above mentioned 
vibrational mode frequency, to any pair of ions. Any superposition of the 2 N states 
of the QC (i.e., the internal states of the ionic system) can then be converted to any 
other superposition by a suitable sequence of one- and two-bit gates. 

In this article we shall analyze the decoherence in the Cirac-Zoller (CZ) QC. A brief 
description of this work has appeared elsewhere [6]. It hardly needs to be said that 
decoherence is a serious limitation to the functioning of any QC. We will be interesed 
in the intrinsic decoherence. Technical difficulties, such as trapping a large enough 
number of ions, proper shaping, phase locking, and timing of the laser pulses, optical 
resolution of individual ions etc., also have the same practical effect as decoherence, 
and may well turn out to be insurmountable by themselves, but that is a separate 
matter. We will mainly discuss the decoherence from the vibrations of the ions. Since 
the discussion in the earlier paper [6] was terse and technical, we will focus in the 
present article on the physical explanation of why ionic vibration is decohering, and 
give a qualitative estimate of the decoherence time. Other discussions specific to the 
CZ scheme are by Plenio and Knight [7], and by Hughes ct al. [8]. General discussions 
of why and how decoherence is detrimental to QC's have been given by Landauer [2] 
and by Unruh [9], among others. 

Intrinsic decoherence in the CZ QC arises from two sources: spontaneous emission 
and ionic vibration. It is easy to estimate the decoherence rate from the first. If the 
spontaneous e) — > \g) decay time for one ion is r s , and we assume that N/2 ions 
are in the excited state on average, we obtain an upper bound of r ra d « 2t s /N for 
the window of time in which any computation must be completed, since a single 
spontaneous decay irretrievably disrupts the wavefunction of the QC as a whole. Our 
estimate ignores effects like superradiance, or changes in the decay of one ion due to 
the proximity of the other ions, but this is justified if the inter ion spacing is larger 
than 2ttc/loq, the wavelength of the g <-> e spectral line. 

A similar estimate of T v ib, the computional time limit imposed by vibrational 
decoherence, is not so easy to obtain. Whatever it is, one simple point should be 
noted now. Since the two mechanisms of decoherence operate independently and in 
parallel, we should add their rates to obtain the total decoherence limit on the useful 
working time of the QC: 

tde C =(r-l+T^y 1 . (1) 

If it should happen that one of the rates, and T~ ih , is much larger than the other, 
this would help us relax the design constraints, as we could then ignore the slow decay 
process to a first approximation. The faster decay may of course still overwhelm us. 

The plan of the paper is as follows. The calculation of T v ib is done in Sec. II. The 
results are discussed in Sec. III. Certain mathematical details are relegated to two 
appendices. 
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II. VIBRATIONAL DECOHERENCE 

The physical origin of vibrational decoherence is as follows. Suppose ion j (sec 
Fig. 1) is not in its equilibrium position. It creates an excess electric field (or electric 
field gradient) on a neighboring ion, i. This excess field alters the evolution in the 
|e), \g) space of ion i from the desired time evolution, and the accumulation of this 
effect causes a decay in the probability that the QC will be in the intended state. 

To qualitatively estimate the decay time, let us denote the longitudinal position of 
the jth ion by Zj , and the deviation from this position by u 3 . Let us further suppose 
that the g <-> e transition is of electric dipole (El) or quadrupolc (E2) type, and 
denote the relevant transition matrix elemenent by d a , where the index a = 1, or 2, 
for the El and E2 cases respectively. Ignoring vector and tensor indices, and denoting 
the ionic charge by q, the change in the eg matrix element of the Hamiltonian for ion 
i is given by 

The key point now is not only that this perturbation is small, i.e., \SVi\ <C 7iwo> 
but that it is also slow. It varies over times set by the periods of the normal modes of 
the ion array, which are much longer than the optical transition time uJq 1 . In other 
words, 

^ ln|^| 



dt 



<w . (3) 



The slowness enables us to treat the perturbation 8Vi adiabatically. Let us map each 
two-state ion onto a spin- 1/2, with e) and \g) being the up and down spin states. 
The ith spin then sees magnetic fields B z = huj and B± — 5Vi (see Fig. 2). Since 
B z » \B±\, the instantaneous precession frequency for this spin is given by 



u4 = {u>* + 5V?/h 2 )V* 



5V 2 

~ W + ^T2^- (4) 

On the other hand, the precession axis for the spin can be taken to be z at all times 
to very good approximation. Thus the time evolution of the spin up and spin down 

states |±) is given by exp (±i f* dt' u>' Qi (t')/2^ |±). A more formal derivation of this 
result is given in Appendix A. 

The phase of the ith ion thus wanders off course by it in a time Tj w n/(uj' 0i — coo). 
Using Eqs. (2) and (4), we can write 

/ / \ 2 \ 

(5) 

The angular brackets above denote some kind of average. The motion of different 
ions is correlated via the normal modes, which we can describe by a density matrix 
at some temperature T. (This is what one means by the statement that the ions 
have been cooled to a temperature T.) Additional correlations due to the systematic 
manipulation of the center of mass mode used to execute the two-bit gates should not 
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be included since this is part of the designed time evolution and does not represent 
a decay. It is therefore completely consistent to have left out this part of the time 
evolution in arriving at Eq. (5). To obtain the best-case answer, we will assume that 
the ion temperature T -C hu> z /kB, and approximate this by T = 0. To obtain an 
order of magnitude, we ignore the details of the normal modes, and simply take 

(ujUk) = -^—Sjk, (6) 

TTiLOt 

where u t is a typical transverse mode frequency and m is the ionic mass. This yields 

q2d * y I . (7 ) 



.,2^2 

rr 1 .- 

2-KhmtjjQUJt , (zi — ZjY 



The rationale for using a transverse mode frequency above is two-fold. First, these 
frequencies are generally higher than the longitudinal ones, and second, it is best to 
choose the states \g) and |e) to be such that the longitudinal modes cannot excite any 
transitions due to a J z selection rule. 

The next step is to obtain the decoherence rate T~ ib for the QC as a whole. One's 
first guess would be that this is obtained by just adding the t~ 1 for all i. The detailed 
calculation of Ref. [6] shows that this is not quite correct. The correct procedure is 
to add the squares and then take the square root. In other words, 

^ b 2 = E-r 2 - (8) 

i 

A qualitative justification for this formula is as follows. The overlap between the 
actual and intended states of the ith spin is better approximated by cos(t/Tj) instead 
of exp(— t/Ti). [See Eq. (21).] The probability P(t) that the QC is in the desired 
state is thus approximately Yli cos 2 (t/Ti) — cxp(— t 2 /r 2 ib ) with r v ib as given above. 
[If we simply added t" 1 , we would obtain an additional factor of iV 1 / 2 in the scaling 
behavior of T~ ih , and Eq. (11) below, e.g., would contain an N instead of the iV 1 / 2 .] 
By combining Eqs. (7) and (8), we obtain a formal answer for the decoherence 
rate. The sums over the lattice positions z t , are at worst, numerical problems. When 
N ^> 1, however, we can evaluate these sums by invoking a continuum approximation 
for the array. This approximation and the lattice sums are discussed in Appendix B. 
The final result can be written as 

vl 2irhmu>ouits 

where so is the minimum spacing between the ions which occurs at the center of the 
array. We can make it apparent that T~ ih is a rate by noting that [sec Eq. (23)], 
q 2 = muj 2 dl, where do is the the trap length scale parameter, and that 

d 2 a Kh/r s k 2 a+1 , (10) 

where fc = w /c. It follows that 

1 N^/doVwl 1 (n) 



T vib T s \S J LUaUJt (k S ) 2a+1 
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Discussion of this result and comparison with the spontaneous decay rate is given 
in the next Section. 

III. DISCUSSION 

To apply Eq. (11) we must take into account that the scaling of tT^ with TV 
depends critically on how the trap operating conditions are varied with TV, and can 
not be naively taken as TV 1 / 2 . If, for example, sq is held fixed as TV is increased, then 
(do/so) 3 ~ N 2 /\nN and tT 1 ~ TV 5 / 2 / In TV. i n this case, however, the longitudinal 
voltage on the trap electrodes, which is proportional to w 2 , varies as (In N)/N 2 . This 
leads to an increase in the total computational time, since the time for executing a 
primitive two-bit gate varies as wj 1 , i.e., as TV. Also, since the longitudinal trapping 
is weaker, the ion array becomes more susceptible to patch voltages on the electrodes, 
and non- linear effects in the trapping potential become more important. If, on the 
other hand, the trap voltages, and therefore, lo z and u t , are held fixed as TV increases, 
then t~-1 ~ jv(8a+i9)/6( m jv)-( 2a + 4 )/ 3 , i. e . ; as jV 9 / 2 (ln TV)~ 2 for an El transition, and 
as TV 35 ' 6 (In TV) -8 / 3 for an E2 transition. In this case, the minimum spacing s varies 
as ~ TV" 2 / 3 , and it may become difficult to resolve the ions optically as is necessary to 
execute the basic gates. Obviously, any intermediate variation is possible by allowing 
both so and uj z to change with TV, and the exact manner in which this is done is thus 
a matter of detailed engineering considerations, which it is premature to discuss. 

To get a numerical estimate of T v ib, we will consider the case of Ba + ions, which 
are particularly favorable from the standpoint of minimizing spontaneous emission 
decoherence. We choose for \g) a state in the ground multiplet 6s 2 Si/ 2 , and for |e) a 
state in the first excited multiplet 5d 2 D 5 / 2 . The frequency ujo = (27r)1.7 x 10 14 Hz. 
Since AL — 2, the 2 D 5 / 2 — > 2 Si/ 2 decay is an E2 process, and the spontaneous decay 
time is r s = 35 s [10]. (There is some uncertainty over this number. Hughes et al. [8] 
take it as 47 s, and it may even be as high as 70 s. We have taken an average from 
Ref. [10].) The index a is 2. Further, we take u z /2n = 100 kHz, and w t /27r = 20 MHz. 
Then do — 14 /im, and for TV = 1000, r vib ~ 10 4 r s , which is surprisingly large. (It 
is even larger in comparison to r ra d = t s /N.) On the other hand, s$ ~ 0.5 /im with 
the same parameters, which runs into the difficulty with optical resolution mentioned 
above. This suggests that a compromise in which lo z is reduced may work better but 
we have not explored this point further. 

That T ra d is so much longer than r v ib for the example chosen above may make 
one wonder if one should not have anticipated this fact. We do not believe so. The 
situation would change completely for larger TV. Indeed for large enough TV, our 
calculation shows that vibrational decoherence will always dominate over spontaneous 
decay decoherence. Secondly, the scaling with TV is quite non-trivial and unexpected. 

We conclude that ionic vibrations are not a significant source of decoherence in the 
original scheme envisaged by CZ, at least for TV < 1000. This should not be taken to 
mean that the problems posed by radiative decoherence by themselves are not serious. 
Indeed, this could well make the whole scheme unworkable. Similarly, the challenges 
of trapping 1000 ions, and of addressing them individually are at the moment quite 
daunting. Nevertheless, our conclusions are encouraging in that they enable us to 
focus on the radiative decay problem. Several authors [7,8] have suggested working 
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with a A system and Raman pulses as a way of dealing with this. An evaluation of 
the vibrational decoherence in this setup remains to be done. 
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APPENDIX A: SPIN-1/2 IN SLOW AND WEAK TRANSVERSE FIELD 

We consider in this Appendix a spin- 1/2 system with the following Hamiltonian 

H=^u a x + f(t)-ff, (12) 

where f(i) = (/ x , f y , 0) has no z component. Further f is small and slow as explained 
in Sec. II. (The field B± introduced in Sec. II is just twice f.) We are interested in 
solving for the time evolution of an arbitrary state for general f(i). 

Denoting the eigenstates of S z with eigenvalues ±1/2 by |±), let us write a general 
state of the spin as 

m)) =«+(t)e-^ ot/2 |+)+7i_(t)e^ t / 2 |-). (13) 
Schrodinger's equation then takes the form 

iu ± = e^Mtfait), (14) 

where f± = f x ±if y . 

Since f± varies very slowly we seek the answer to Eq. (14) in the form 

u±(i) = a±(t) + f3±(t), (15) 

where a± and (3± are fast and slow parts, the latter varying little over a period 2tt/ui , 
and the former averaging to zero over several such periods. 

Substituting Eq. (15) in (14), and separating the fast and slow parts, we obtain 

zd ± = e ±iuot Mt)^(t), (16) 
± = e^f^a^t). (17) 

To integrate Eq. (16), it is a good approximation to treat the slowly varying functions 
f± and (3± as constants. In this way, we obtain 

a± = Tu>o 1 e ±iaJot Mt)(3 T (t). (18) 

We now put this solution in Eq. (17), and average the resulting equation over several 
periods 27r/u - This yields 

i/ j ± = ±MJ! / a t ( t ). (19) 

LJq 
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Integrating this, we obtain (3±(t) — exp(=Fz$(t))/3±(0), where 

$(*) = f dt'lffl!. (20) 
Jo 

We can take u± « /3±, since a± <C |/3±|. The quantity <&(t) is precisely the excess 
angle through which the spin precesses about the z axis as discussed in Sec. II, i.e., 
it is half the difference between J Q u>' oi (t') dt' and uiot. 

Suppose the initial state of the spin is 2~ 1 / 2 (|+) + |— )), i.e., u±(0) = 2~ 1//2 . Let 
us denote the state at time t that would be obtained in the absence of the transverse 
field f by \ipo(t)). In the context of the QC, the states \"4>o(t)) and \if>(t)) are analogous 
to the states of the ideal and actual QC, without and with decoherence, respectively. 
The extent of decoherence is given by the overlap 

<^o(t)|^(t))=cos($(t)). (21) 

A more careful calculation shows that $(i) also has a Berry phase part z • (f xf). 
For the problem of interest to us, this is much smaller then the dynamical phase and 
may be neglected. Further, the adiabatic approximation breaks down due to secular 
effects for t ~ <x>o/|f| 2 . It is not hard to see that this breakdown time is much larger 
than T v ib, and so the approximations (4) or (20) are completely adequate. 



APPENDIX B: CONTINUUM APPROXIMATION FOR ION ARRAY 

We wish in this Appendix to quantitatively understand the structure of the linear 
array of trapped ions when N 1. We do this via a continuum approximation based 
on the expectation that the local spacing s{zi) between ions in the vicinity of ion i 
will vary slowly with i. Our goal is to find the function s(z), where we regard z as 
a continuous variable. We will also find how the total length of the array, 2L, varies 
with N. 

A simple-minded argument for s(z) is as follows. Consider the Coulomb forces 
on an ion at position z from its nearest neighbours to the left and right, which 
we take to be at positions z — s_ and z + s + respectively. The net force is then 
q 2 (sZ 2 — s+ 2 ) w 2q 2 s~ 2 (ds/dz), where we have approximated s + — s_ by s(z)(ds/dz). 
The force from successively distant pairs of neighbours is smaller than this expression 
by factors of 4, 9, 16, etc., since the distances are approximately doubled, tripled 
and so on. Thus the net Coulomb force is (n 2 q 2 /3s 2 )(ds/dz), since ^2 n n~ 2 = ir 2 /6. 
Equating this to the opposing spring force muj 2 z from the trapping potential, we 
obtain 

(22) 



3s 2 (z)dz dg' 
where 

d Q = (q'/mLo 2 ) 1 / 3 (23) 

is a natural length scale for the trap. (It is easy to show that the ion spacing is of the 
order of do for 2 or 3 ions in the trap.) Denoting the total length of the array by 2L 
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and placing the center at z = 0, integration of Eq. (22) gives 



wrwr-^ iL '-^ (24) 

By balancing the forces on the ion at the end of the array in the same way as was 
done above, we obtain s(L) w ■jr(d 3 t /6L) 1 ^ 2 . We can thus ignore s _1 (£) compared to 
?>L 2 /2iT 2 dl in Eq. (24), which yields 

s(z) = s (l-z 2 /L 2 )-\ (25) 

where so = s(0) = 2-K 2 d 3 ) /L 2 is the minimum ion spacing (attained at z = 0). 

To obtain an expression for L(N), let us denote the ion number at position z by 
n(z). Then, integration of the approximate relation dn/dz — l/s(z) using Eq. (25) 
gives L = d ('ir 2 N/2) 1 / 3 . This implies so ~ N~ 2 I 3 . The mean spacing can be found 
to vary as N~ 2 / 3 In TV. 

It is clear that the above argument does not treat the ends of the chain properly, 
and also underestimates the Coulomb forces due to the more distant neighbors. A 
more sophisticated approach is due to Dubin [11], who treats the ion array as a fluid 
of total charge Nq. In a harmonic trapping potential, the solution to this problem is 
known: the fluid forms a blob of uniform charge density in the shape of an ellipsoid 
of revolution of total volume AirNd^ provided ui t 3> u> z . If the semi major axis of 
this ellipsoid is L, the semi minor axis is therefore (ZNdl/L) 1 / 2 . The inverse spacing 
l/s(z) is clearly the number of charges per unit length along the major axis of the 
ellipsoid, and this in turn is given by its cross-sectional area. In this way we obtain 



a(z) 4 L V L 2 

This is identical to our approximate form (25), but s is different. 

We still need to find the length L. The fluid approximation breaks down over 
here, as the answers it gives depend on the ratio u) z /u t , which is clearly wrong as 
long as the linear structure is stable. Dubin uses a local density functional theory to 
estimate the correction to the Coulomb energy due to the discreteness of the array, 
and minimizes the sum of this correction, the fluid drop self energy, and the trapping 
potential energy, with respect to L. The result for L is then independent of uj z /u} t 
and is given by 

L 3 = 3-/Vln(0.87V)c?Q. (27) 

(The 0.8 is actually 6e 7 " 13 / 5 , with 7 being Eulcr's constant.) This result differs from 
our previous one by logarithmic factors. We also obtain 

s = 4X/37V = 1.92iV- 2 / 3 [ln(0.8iV)] 1/3 do. (28) 

This result should also be compared to that of Hughes et al. [8], who find on the 
basis of a numerical fit that so = 2N~°- 56 do. Since for moderate values of N, the 
logarithmic factor in Eq. (28) will have the effect of increasing the apparent exponent 
of N, these two results are quite comparable. We do not know if Hughes ct al. did 
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numerics for large enough N to discern the presence or absence of logarithmic factors, 
so it is hard to say which is better. 

The above results can be used to perform the sums over lattice positions that 
appear in Eqs. (7) and (8). There are two types of sums. The first, 



can be very simply evaluated as 



5 n (<)«25»(^)f;4=#l. (30) 

^J™ S n (Zi) 

This approximation is actually fairly good for all i except very close to the ends, since 
the exponent n is at least 3 or 4 in all cases that we encounter. 

The second type of sum is T n — J2i s~ n (zi). This can be approximated by an 
integral. Writing Ai w dz/s(z), we obtain 

L 

Tn= ^^{z7)~ I ^{z)- (31) 

1 -L 

With s(z) given by Eq. (25), the integral is elementary, and we obtain 

^M^T' (32) 

where the last form comes from an asymptotic formula for p(n + 2, 1/2). 
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LIST OF FIGURES 

Fig. 1. Schematic of the array of trapped ions in the Cirac-Zoller quantum computer. 

Fig. 2. Equivalent magnetic fields acting on the internal states of the ith ion in the 
spin representation. 
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